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I. INTRODUCTION 

In a recent paper [1] equivalence between isokinetic and 
frictionless thermostats has been discussed heuristically, 
leaving aside several difficulties on the understanding of 
the classical dynamics of systems of infinitely many par- 
ticles. Such an understanding is, however, a necessary 
prerequisite because strict equivalence can be expected 
to hold only in the thermodynamic limit. In this paper 
we proceed along the same lines, comparing the isoen- 
ergetic and the frictionless thermostats, and study the 
conjectures corresponding to the ones formulated in [1] 
for isokinetic thermostats, obtaining a complete proof of 
equivalence in 1, 2, 3-dimensional systems with various 
geometries. 

In Sec. II the class of models to which our main re- 
sult applies is described in detail. The main result is 
informally quoted at the end of the introduction after 
discussing the physics of the models; a precise statement 
will be theorem 1 in Sec. IV and it will rely on a property 
that we shall call local dynamics: the proof is achieved by 
showing that in the models considered the property holds 



as a consequence of the theorems 2-10, each of which is 
interesting on its own right, discussed in the sections fol- 
lowing Sec. IV and in the appendices. 

A classical model for nonequilibrium statistical me- 
chanics, e.g. see [2], is a test system in a container fig, 
and one or more containers Qj adjacent to it and enclos- 
ing the interaction systems. 

A geometry that will be considered in dimension d = 
2, to fix the ideas, can be imagined (see Fig.l, keeping 
in mind that it is just an example for convenience of 
exposition and which could be widely changed) as in the 
first picture in Fig.l; the second picture illustrates the 
only gemoetry that we shall consider in dimension d = 3: 




• • : Xy, Xy) 

Fig.l: The 1 + u finite boxes fij n A, j = 0, . . . , z^, are marked 
Co , Ci , . . . , Ci/ and contain No, Ni, . . . , N^, particles out of the 
infinitely many particles with positions and velocities denoted 
Xo, Xi, . . . , X^, and Xo, Xi, . . . , X^, respectively, contained 
in j > 0. The second figure illustrates the special geometry 
considered for d = 3; here two thermostats, symbolized by the 
shaded regions, , occupy half-spaces adjacent to Qo ■ 

Referring, for instance, to the first of Fig.l: 

(1) The test system consists of particles enclosed in a 
sphere ilo — S(I?o) of radius Do centered at the origin. 

(2) The interaction systems consist of particles enclosed 
in regions which are disjoint sectors in m."^. I.e. disjoint 
semiinfinite "spherically truncated" cones adjacent to fio, 
of opening angle ujj and axis k^: VLj = e K'', |C| > 
Dq,\- kj > Id coswj}, j ^l,...,v. 

The initial configurations x of positions and veloci- 
ties will be supposed to contain finitely many particles 
in each unit cube. Thus the test system will consist of 
finitely many particles, while the interaction systems are 
infinitely extended. 

We shall suppose that the forces acting on the particles 
are due to a repulsive pair interaction of radius and 
to external repulsive interactions acting within a distance 
from the boundaries of the containers plus in the not 
frictionless models some thermostatting forces (of Gaus- 
sian type), [3]. Furthermore on the test system may act 
a nonconservative "stirring" force ^. 

The motion starting from x must be defined by first 
regularizing the equations of motion (which are infinitely 
many and therefore a "solution" has to be shown to exist) 
"approximating" them with evolution equations involv- 
ing finitely many particles. 

There is wide arbitrariness in the choice of the reg- 
ularization: and the ambiguity should disappear upon 
regularization removal. 
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For instance a first regularization, that we call elas- 
tic reqularizafAon, could be that only the (finitely many) 
particles of the initial data x inside an artificial finite 
ball A = E(r) of radius r > Dq will be supposed moving 
and will be kept inside A by an elastic reflection bound- 
ary condition at the boundary of A. The particles of x 
located outside the container rioUUj>o(i^j H A) are imag- 
ined immobile in the initial positions and influence the 
moving particles only through the force that the ones of 
them close enough to the boundary of A exercise on the 
particles inside A. 

A second alternative regularization, that we shall call 
open regularization, is obtained by letting only the par- 
ticles initially inside A move while their motion is not 
influenced by the particles external to A and they are 
even allowed to exit the region A. 

In the ^thermodynamic limit " , which will be of central 
interest here, the ball A grows to oo and the particles that 
eventually become internal to A start moving: in other 
words we approximate the infinite volume dynamics with 
a finite volume one, called A-regularized, and then take 
an infinite volume limit. 

In the A-regularized evolution the energy in the re- 
gion Q,j n A, J > 0, will in gcnciral change in time. The 
isoenergetic thermostat is defined by adding "frictional" 
forces —ajQi on all particles {qi,qi) in flj fl A where aj 
is chosen in such a way that the total energy in n A 
is constant in time, see the next section for details. The 
thermostatted evolution is then the evolution when such 
frictional forces are added in each ftj C] A, j > 0. 

Isoenergetic models of the kind considered here have 
been studied also in simulations aiming at checking the 
thermostats "efficiency", i.e. the possibility of a bound- 
ary heat exchange sufficient to allow reaching stationarity 
in systems with many particles, [4]. 

The essential physical requirement that the thermo- 
stats should have a well defined temperature and density 
will be satisfied by an appropriate selection of the ini- 
tial conditions. The guiding idea is that the thermostats 
should be so large that the energy that the test system 
transfers to them, per unit time in the form of work Qj , 
is acquired without changing, at least not in the ther- 
modynamic limit, the average values of the densities and 
kinetic energies (i.e. temperatures) of the thermostats in 
any finite observation time > 0. 

To impose, at least at time and in the thermodynamic 
limit, the requirement the initial data will be sampled 
with probability /io where is the product of extremal 
DLR distributions in each ftj with distinct temperatures 
and chemical potentials, see the next section for details. 

Main result: In the thermodynamic limit, the ther- 
mostatted evolution, within any prefixed time interval 
[0,Q], becomes identical to the frictionless evolution at 
least on a set of configurations which have iiq -probability 
1 with respect to the initial distribution no, in spite of the 
non stationarity of the latter. In the same limit also the 
frictionless evolution with open or elastic regularization 



become identical. 

This is proved after theorem 6 in Sec. VII. 

The present paper relies heavily on results available 

in the literature, [5 8], but the results cannot be directly 
quoted because (minor) modifications of them are needed 
for our purposes. Therefore this paper is self-contained 
and the results that could be also found by a careful study 
of the literature are proved again from scratch: however 
we have confined such derivations in the Appendices A-F 
with the exception (for the sake of clarity ot the argu- 
ments) of Sec. V, VI; the Appendices G-L give technical 
details about the new methods, explained in Sec.VII and 
in the Appendices G-J. 

The strong restriction on the thermostats geometry in 
(i = 3 is commented in Appendix D after the proof of 
lemma 2, where it is used. 



II. THERMOSTATS 

A configuration x is thus imagined to consist of a 
configuration (Xo,Xo) with Xq contained in the sphere 
T,{Dq), delimiting the container Qq of the test systems, 
and by configurations (Xj,Xj) with Xj C Ctj, Xj S 
R''),j = l,...,u: 

Phase space: Phase space % is the collection of locally 
finite particle configurations x = {. . . ,qi,qi, . . 

X = (Xo, Xo, Xi, Xi, . . . , X,, X,) =^ (X, X) (2.1) 

with Xj c Oj, qi e IR"^; in every ball S(r') of radius r' 
centered at O, fall a finite number of points ofX.. 

The particles of x located outside A will be regarded 
as immobile or even as non-existing (depending on the 
choice of elastic or open regularization considered). It 
will be convenient to suppose that the regularization re- 
gion A is a ball A„ of radius 2"r(^, n = no, no -I- 1, . . 
with no large so that 2"''r,^ > Do-\-r^ so that A„ contain 
the test system and the particles interacting with it. 

The particles are supposed to interact with each other, 
via a potential ip, and with the non artificial walls {i.e 
those of the containers ilj), via a potential ip: 

Interaction: Interparticle interaction will be through a 
pair potential with finite range r^ smooth, decreasing 

and positive at the origin. The walls of the containers Qj 
are represented by a smooth decreasing potential V' ^ o/ 
range r^p <C r^p and diverging as an inverse power of the 
distance to the walls UjdO,j, while the (artificial bound- 
ary) OA will be imagined as a perfectly reflecting elastic 
barrier in the case of elastic boundary conditions or as a 
boundary perfectly transparent to the moving particles. 

Hence the potential (f is superstable in the sense of [9] : 
a property that will play an important role in the fol- 
lowing. The value of the potential ip at midrange will be 
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denoted and < (p'^= <f{^) < "^(0); the wah potential 
at distance r from a wall will be supposed given by 



where f A^^ As Tj(x) =^ /C^- av-j; - ^^^j "J 
chosen so that Kj^A(X.j) + J7j,A(Xj) = Ej^\ are exact 
constants of motion: the subscript A will be omitted un- 
less really necessary. A more general model to which the 
analysis that follows also applies is in [10] . 

(5) The forces #(Xo) arc, positional, nonconservative, 
smooth "stirring forces" (possibly absent). 

(6) In the case of A-regularized thermostatted dynamics 
we shall consider only initial data x for which the kinetic 
energies Kj^^i^j) of the particles in the f2j n A's are 
> for all large enough A. Then the time evolution is 
well defined for t < tA{x) where t^ix) is the maximum 
time before which the kinetic energies remain positive 
(hence the equations of motion remain defined for t < 
tA{x) because the denominators in the aj stay > 0; see 
Appendix I for technical details). 

The equations with o = 1 will be called A-regularized 
isoenergetically thermostatted because the energies Ej = 
Kj + Uj stay exactly constant for j > and equal to their 
initial values Ej. The equations with a = in Eq.(2.3) 
will be considered together with the above and called the 
A-regularized frictionless equations. 

Remark that Qj is the work done, per unit time, by 
the test system on the particles in the j-th thermostat: 
it will therefore be interpreted as heat ceded to the j-th 
thermostat. 

It will be convenient to consider also the open bound- 
ary condition at least in the frictionless thermostats case: 
the equations of motion are immediately written. 

To impose, at least at time and in the thermody- 
namic limit, the requirement that the thermostats should 
have a well defined temperature and density the values 
Nj,Ej,j > 0, will be such that |^^^^| ^^^> Sj and 

,„ . I ^e,-: with , e,- > fixed in a sense that is 

|UjnA| A— 7-co J J 

specified by a choice of the initial data that will be stud- 
ied, and whose physical meaning is that of imposing the 
values of density and temperature in the thermostats, for 
J >0. 

Initial data: The probability distribution fiQ for the ran- 
dom choice of initial data will be, ifdx'^= 0^=0 '~~nT^' 
the limit as Aq —>■ oo of the finite volume grand canonical 
distributions on H 



(Xj) 1i/'fX2 and aj are 



. r^ 
r < 



(2.2) 



with a > and r equal to the distance of q to the wall; 
for larger r it continues, smoothly decreasing, reaching 
the value at r = r^. The choice of ^f) as proportional 
to ipo limits the number of dimensional parameters, but 
it could be made general. The restriction <Si r^ is not 
necessary: it has the physical interpretation of making 
easier the interaction between particles in fig and parti- 
cles in Uj>of2j and, therefore, transfer of energy between 
test system and thermostats. 

Particles in ilo interact with all the others but the par- 
ticles in Vlj interact only with the ones in J7j U i7o: the 
test system in interacts with all thermostats but each 
thermostat interacts only with the system, see Fig.l. 

The A-regularized elastic boundary condition equa- 
tions of motions (see Fig.l), aside from the reflecting 
boundary condition on the artificial boundary of A, con- 
cern only the particles in JIq U Uj>o(rij- n A) and will be 



mXoi = - 5iC/o(Xo) - 5.t^o,j (Xo, X,) + ^^(Xo) 
mXji = - diUj{y.j) - aiC/o,j(Xo, X,) - a ajXji (2.3) 



(1) where the parameter a will be a = 1 or a = depend- 
ing on the model considered; 

(2) the potential energies Uj{Xj), j > and, respec- 
tively, C/o,j(Xo,Xj) denote the internal energies of the 
various systems and the potential energy of interac- 
tion between the system and the thermostats; hence for 
Xj c n A the Uj's are: 



[/,(X,)= ^ i;{q)+ ^ ip{q-q') 
C/o,j(Xo,Xj) = J2 fil-l'h 

9GXo,9'eXj 



(2.4) 



(3) the first label in Eq.(2.3), j — or j ~ 1, . . . ,1/, re- 
spectively, refers to the test system or to a thermostat, 
while the second labels the points in the corresponding 
container. Hence the labels i in the subscripts {j, i) have 
Nj values and each i corresponds (to simplify the nota- 
tions) to d components; 

(4) The multipliers aj are, ioi j = 1, . . . ,u, 



def 

a-j = 



Qj 



dNjkBTj{x)/m' 



with 



def 



(2.5) 



Xj ■ djUo,j{XQ,Xj), 



tioAoidx) = const e--^o-^o(^) dx, with (2.6) 

Ho,Ao{x) =^ 5];3,(i^,-,Ao(x) - A,%Ao + Uj,Ao{x)) 

i=o 

Remarks: (a) The values /3o = j^j^ > 0, Ao € M, are also 
fixed, although they bear no particular physical meaning 



4 



III: Notations and sizes 



because the test system is kept finite. 

(b) Here A = (Ao,Ai,...A^) and T = (Tq, Ti, . . . , T^) are 
fixed chemical potentials and temperatures, and Aq is a 
ball centered at the origin and of radius Dq. 

(c) The distribution fiQ is a Gibbs distribution obtained 
by taking the "thermodynamic limit" Aq — >■ oo. The the- 
ory of the thermodynamic limit implies the existence of 
the hmit distribution /xq, either at low density and high 
temperature or on subsequences, [11]. In the second case 
(occmring when there arc phase transitions at the chosen 
values of the thermostats parameters) boundary condi- 
tions have to be imposed that imply that the thermostats 
arc in a pure phase: for simplicity such cxcciptional cases 
will not be considered; this will be referred to as a "no- 
phase transitions" restriction. 

(d) Ao should not be confused with the regularization 
sphere A: it is introduced here and made, right away, oo 
only to define /io- 

(e) Notice that /io is a product of independent equilibrium 
Gibbs distributions because Hq does not contain the in- 
teraction potentials Uqj. 

(f) The proofs extend to /Lto's which are product of DLR 
distributions in each container (which are not necessarily 
extremal). 



III. NOTATIONS AND SIZES 



Initial data will be naturally chosen at random with 
respect to /io- Let the "pressure" in the j-th thermostat 

be defined by pj A; Ao) =^ ^3 |a/nAo| ^jM iP^ ^) with 



00 

E 

JV=0 



-l3{-\N+Kj{xN) + Uj{xM)) 



(3.1) 



where the integration is over positions and velocities of 
the particles in Ao n ^Ij. Defining p{(3,\) as the ther- 
modynamic limit, Ao 00, of pj(/3. A; Ao) we shall say 
that the thermostats have densities 5j, temperatures Tj, 
energy densities Cj and potential energy densities Uj , for 
j > 0, given by equilibrium thermodynamics, i.e.: 



if X — )• SfX denotes the solution of the equations of motion 

(for the considered model). The variation (t(.x), per unit 
time of a volume element, in the sector of phase space 
containing Nj > particles in fljOA, j = 0,1, . . . ,1/, can 
be computed and is, under the A„-regularized dynamics 
and for the elastic boundary conditions models. 



a(a;)-^/3,g„ a = (3.4) 

j>0 



j>0 



def 

as it follows by adding the time derivative PoQo = 
^o{Ko + Uq) to the divergence of Eq.(2.3) (regarded as a 
first order equation for the g's and g's) using the expres- 
sion in Eq.(2.5) for aj. 

Remarks: (1) The expressions for (t{x) will play a central 
role in our approach and one can say that the key idea 
of this work is to control the size of |cr| and through it 
show that "singular" events like all particles come close 
to a simultaneous stop or accumulate to a high density 
somewhere or some of them acquire very high speed or 
approach too closely the walls (hence acquiring huge po- 
tential energy) will have zero /io-probability. Such events 
have zero /to-probability at time and the variations with 
time of Ho are controlled by |a|. Naturally, in view of 
Eq.(3.4), we shall call such estimates "entropy bounds. 

(2) The relation /3o (Ko + Uo) = /3o (* • Xq - E,>o (Uoj - 
Qj)) is useful in studying Onsager reciprocity and Green- 

def 

Kubo formulae, [12]. Notice that Kq + Uq = Qo is also 
- Ej>o Xo • axo?7o,,(Xo, X,) + # • Xo. 

(3) It is also interesting to consider isokinetic ther- 
mostats: the multipliers aj are then so defined that 
Kj{x) is an exact constant of motion: calling its value 
^NjkBTj{x) the multiplier aj becomes 



aj{x) 



def 



d NjkBTj{x)/m' 



(3.5) 



(3.2) 



_ dl3jp{l3j,Xj) 
^3 - ^j^j, 



which are the relations linking density Sj, temperature 
Tj = {kBi3j)~^, energy density and potential energy 
density Uj in a grand canonical ensemble. 

In general the A regularized time evolution changes 
the /xo-measure of a volume element in phase space by 
an amount related to (but different from) the variation 
of the Liouville volume because, in general, 



d Ho{Stdx) I 



dt 



IJ-oidx) 



(3.3) 



with Qj defined as in Eq.(2.5). They have been studied 
heuristically in [1]. 

The open boundary conditions will only be considered 

for the A„-regularized frictionless dynamics (A„ is de- 
fined in the paragraph following Eq.(2.1)) and denoted 

X — >■ 5|"'°'x: in this case only the particles initially in- 
side A„ will move freely allowed to cross 9A„ unaffected 
by the particles initially outside A„ which will keep their 
positions fixed. 

The A„-regularized motions with elastic boundary 
conditions will be denoted x si^'^'^x, a — 0,1. 



Let S^"'^ X or X be the infinite volume dynamics 



limn 



0,1, or lim„- 



tt(".0) 



X provided 



the limits exist. We shall often use the notations 
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(3.6) 



Remarks: (1) In the frictionless case the existence of a 
solution to the equations of motion poses a problem only 
if we wish to study the A„ oo limit, i.e. in the case in 

which the thermostats are infinite: for A„ finite S'["'°''x 
is well defined. While for tlie clastic reflections at the 
artificial regularization boundary 9A„ it is shown in [7] 
that the dynamics is also well defined with /xo-probability 
1. 

(2) In the thermostatted case the kinetic energy appear- 
ing in the denominator of aj, see Eq.(2.5), can be sup- 
posed to be > with /Uo-probability 1 at f = 0. However 
it can become later at tf^^{x) (see item (6), p. 3, and 
the example at end of Sec. IV). In the course of the anal- 
ysis it will be proved that with /xq -probability 1 it is 
tA„ (a;) „^oo> oo; therefore S^^'^^x is eventually well de- 
fined. 

(3) It will be shown that the Umits {t) and a;(°) [t) ex- 
ist and are identical (as expected). It will become clear 
why it is necessary to consider first the frictionless ther- 
mostats with open boundary conditions. 

We shall denote {S'f''°'^x)j or a;^"'"^(t) the positions 

and velocities of the particles of S'j"'"''.t in a corre- 
sponding notation will be used for positions and velocities 



of S. 



{n,0) 



X . 



lixfi'^'it) denote the pairs of positions and velocities 

((z|"'°''(i): ?■"'"•'(*)) with qi G Then a particle with 
coordinates {qi,qi) at t = in, say, the j-th container 
evolves between collisions with the regularization bound- 
ary A„ (if any), see Eq.(2.3), as 



Qiit) =%(0)+ / qiit')dt' 
Jo 

q,{t) =e-/o»"^(*')'^*'g,(0) 

+ /* e" ""^ /i(a;("'''H*")) dt" 
Jo 



(3.7) 



where m/, = -9,, (C/, (X, (i)) + C/o,i(Xo(t), X, (t))) 



SjoM^oit)) and Xj{t) denotes xj"'"^(i) or X.p"'{t). 
Here Uj is defined as in Eq.(2.4); ao = 0. 

The open boimdary evolution is described by Eq.(3.7) 
with a = and with the appropriate interpretation of fi. 

The first difficulty with infinite dynamics is to show 
that the number of particles, and their speeds, in a finite 
region A remains finite and bounded only in terms of the 
region diameter r (and of the initial data): for all times 
or, at least, for any prefixed time interval. 

We shall work with dimensionless quantities: therefore 
suitable choices of the units will be made. If 9 is prefixed 
as the maximum time that will be considered, then 



(fo : (energy scale), r<p ; (length scale) 
Q : (time scale), vi - 



2ip(0) (3-8) 
(velocity scale) 



m 

are natural units for encirgy, Icingth, time, velocity. 

It will be necessary to estimate quantitatively the size 
of various kinds of energies of the particles, of a config- 
uration X, which are localized in a region A. Therefore 

introduce e{q, q) + V-'(9)) / ^o and, for any region 

A, the following dimensionless quantities: 



(a) Na{x), Nj^A{x) the number of particles of x 
located in A/Oq or, respectively, A n Qj 



(3.9) 



def 

(6) e^ix) = max e{qi,qi) 

(c) Ua{x) = ^ ^ ifi{qi - qj)/^o 

{d) Va{x) = max \qi\/vi 

g,eA/no 



The symbol B{£,,R) will denote the ball centered at 
^ and with radius i?r<^. With the above notations the 
local dim,ensionless energy of the thermostat particles in 
B{£,,R) will be defined as 



W{x;^,R)''^^ ^ 



E ( 



rnq^ 



ip(0) ^ \ 2 



1 



(3.10) 



Let log+2:''= max{l,log2 \z\}, gQ{z) = (log+z)'^ and 



C(^[x) = sup 



sup 



Wix;^,R) 
Rd- 



(3.11) 



If % is the space of the locally finite configurations {i.e. 

containing finitely many particles in any finite region), let 
C H he the configurations with 

W'd'Xoo, (2)^>lg (3.12) 

for all A = B{0,L) large enough and for Sj,Tj, given 
by Eq.(3.2); here Nj,\, Uj,\,Kj,\ denote the number of 
particles and their potential or kinetic energy in ilj D A. 

Remark: Notice that the lower bound in (2) of Eq.(3.12) 
is half the value of the average kinetic energy in the initial 
data and it will become clear that any prefixed fraction 
of the average kinetic energy could replace 1/2. Each set 
Hi; has /XQ-probability 1 for ^ > 1/rf, see Appendix A. 
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IV. EQUIVALENCE: ISOENERGETIC VERSUS 
FRICTIONLESS 

Adapting a conjecture, in [1] , we are led to expect that 
the motions considered satisfy the property: 

Local dynamics Given 9 > 0, for t G [0, ©] and with 

fiQ-probability 1: 

(1) The limits x^''\tf= lim x^''^'-'\t) , x^^^t)"^^ lim 

n— >-oo n— >-oo 

2;(">o)|'^-j ^"thermodynamic limits") exist for all t < Q 
and a = 0, 1. 

(2) Fort< Q, a;("'iH*) saMsfi.es the second of Eq.(3.12). 

(3) The functions t xJ"'>{t) and t x'^"\t) solve 
uniquely the frictionless equations in a subspace of H to 
which also x^^\t) belongs (explicit, sufficient, bounds are 
described in theorem 5). 

Remarks: (a) The limits of as A„ 

oo, arc understood in the sense that for any ball A whose 
boundary does not contain a particle of x''^\t) the la- 
bels of the particles of x'''^\t),x^^\t) and those of the 
particles in x^"'"^ (i), a;^"'°^ (i) which are in A are the 
same, and for each i the limits of ((/f"'"'' (i), e[[^''^\t)) and 

exist and are continuous, together 
with their first two derivatives. 

(b) Uniqueness in item (3) can be given several meanings. 

The simplest is to require uniqueness in the spaces 
for Q>lld fixed: and theorem 11, in Appendix F, shows 
that for = 1,2 one could suppose such simpler property. 
However our statement is more general and we have left 
deliberately undetermined which subspace is meant in 
item (3) so that the determination of the subspace has to 
be considered part of the problem of establishing a local 
dynamics property. The generality becomes relevant in 
studying the case c? = 3, where even in equilibrium a 
proof that the evolution of data in T-Lq remains in the 
same space is lacking. The local dynamics property in 
d = 3 is implied by theorem 9, in Appendix E. 

(c) Recalling the characteristic velocity scale (namely 
v\ = y'2(p(0)/m), the initial speed of a particle located 

in 5 e M'', is bounded by f i y^f i /d9i /dil/TipY^"^', and 
the distance to the walls of the particle located at q is 
bounded by (Si/d9i/d{<llri,Yy^'°' r^. 

Hence for large they are, respectively, bounded pro- 
portionally to [(log |(7|/r;p)3] 2 and [(log Igl/rip)^]^- : this 
says that locally the particles have, initially, a finite den- 
sity and reasonable energies and velocity distributions 
(if measured on boxes of a "logarithmic scale"). The 
theorem 11 in Appendix F will show that this property 
remains true for all times, with /xq probability 1. 

(d) An implication is that Eq.(3.7) has a meaning at time 
t = Q with /XQ-probability 1 on the choice of the initial 
data x, because £{x) < oo. 

(e) The further property that the thermostats are effi- 
cient: i.e. the work performed by the external non con- 



servative forces is actually absorbed by the thermostats 
in the form of heat Qj , so that the system can eventually 
reach a stationary state, will not be needed because in a 
finite time the external forces can only perform a finite 
work (if the dynamics is local). 

(f) It should also be expected that, with /io^probability 1, 
the limits as A — >• oo of item (2), Eq.(3.12), should exist 
and be equal to ^ for almost all t > respectively: 

this is a question left open (and it is not needed for our 

purposes). 

Assuming the local dynamics property, equivalence, 
hence the property x^^^ {t) = x^^'> (t) = a;^^^ (t) for all finite 
t, can be established as in [1]. This is recalled checking, 
for instance, x^^\t) = x^^\t), in the next few lines. 

In the thermostatted A„-regularized case it is 

|aj(a;)| = ^ — ^ (4.1) 

The force between pairs located in , Qj is bounded by 

def 

F = max |9(p(g)|; the numerator of Eq.(4.1) can then be 
bounded by FNq \/N&^ 2Kj /m where A^o is the number 

of particles in Cq and No bounds the number of thermo- 
stat particles that can be inside the shell of radii Do, 
-Do + Tip for < t < (by Schwartz' inequality). 

Remark that the bound on Nq exists by the local evo- 
lution hypothesis (see (1) and remark (a)) but, of course, 
is not uniform in the initial data x. Hence for < t < Q 
and large enough n Eq.(3.12) yields: 

" V2iir,-,A„(x(i.A„)(t)) - ^|r!,nA„|(5,(i/2/3/ 

(4.2) 

Letting A„ — > oo it follows that aj ^ _^^> 0. 

Taking the limit A„ ^ oo of Eq.(3.7) at fixed i, this 
means that, with iiq -probability 1, the limit motion as 
A„ — ^ oo (with l3j,Xj, j > 0, constant) satisfies 

qi{t) = qi+ [ qi{t')dt', qi{t) = qt + [ fi{t")dt" (4.3) 
Jq Jo 

i.e. the frictionless equations; and the solution to such 

equations is uniqiie with probability 1 (having again used 
assumption (3) of the local dynamics). In conclusion 

Theorem 1: // the dynamics is local in the above sense 
then in the thermodynamic limit the thermostatted evolu- 
tion, within any prefixed time interval [0, O], becomes the 
frictionless evolution at least on a set of configurations 
which have probability 1 with respect to the initial distri- 
bution /io ; in spite of the non stationarity of the latter. 

Suppose, in other words, that the initial data are sam- 
pled with the Gibbs distributions for the thermostat par- 
ticles (with given temperatures and densities) and with 
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an arbitrary distribution for the finite system in ilo with 
density with respect to the Liouvillc volume (for instance 
with a Gibbs distribution at temperature (fcs/Jo)"^ and 
chemical potential Aq as in Eq.(2.6)). 

Then, in the thermodynamic limit, the time evolution 
is the same that would be obtained, in the same limit 
A„ — 00, via a isoenergetic thermostat acting in each 
container flj fl A„ and keeping its total energy (in the 
sector with Nj particles) constant and with a density 
equal (asymptotically as A„ — >■ oo) to e^. 

The difficulty of proving the locality property (2) can- 
not be underestimated, although it might seem, at first 
sight, "physically obvious" : the danger is that evolution 
implies that the thermostat particles grind to a stop in 
a finite time converting the kinetic energy entirely into 
potential energy. The consequence would be that aj be- 
comes infinite and the equations ill defined. 

As a consequence it is natural to expect, as stated 
in the local dynamics assumption, only a result in /io~ 
probability. This can be better appreciated considering 
the following counterexample, in the frictionless case. 

Consider an initial configuration in which particles are 
on a square lattice (adapted to the geometry): regard the 
lattice as a set of adjacent tiles with no common points. 
Imagine that the particles at the four corners of each tile 
have velocities of equal magnitude pointing at the center 
of the tile. Suppose that the tiles sides are > r^. If (fi{0) 
is large enough all particles come to a stop in the same fi- 
nite time and at that moment all kinetic energy has been 
converted into potential energy: at time all energy is 
kinetic and later all of it is potential. Certainly this ex- 
ample, which concerns a single event that has, therefore, 
probability in jiQ, shows that some refined analysis is 
necessary: the thermostatted evolution a;("'^^(t) might 
be not even well defined because the denominator in the 
definition of aj might become 0. 

It should be stressed that the thermostats models con- 
sidered here preserve, even at finite A„, an important 
symmetry of nature: time reversal: this certainly ex- 
plains the favor that they have received in recent years 
in the simulations. 

A corollary will be that neither the frictionless motion 
nor the dissipative thermostatted motions conserve phase 
space volume (measured with no), but in both cases the 
entropy production rate coincides with the phase space 
volume (measured by //q) contraction and, at the ther- 
modynamic limit, is identical in the two cases. 

In Appendix G we have discussed the very simple case 
in which there is no direct mutual interaction between 
the thermostats particles, free thermostats, and the walls 
are reflecting. Reading this simple case may help under- 
standing the organization of the analysis in the interact- 
ing cases. 



V. ENERGY BOUND FOR FRICTIONLESS 
DYNAMICS 

Proof of the local dynamics property will require con- 
trolling the maximal particles speeds, the number of par- 
ticles interacting with any given one as well as their num- 
ber in any finite region. This will be achieved by proving 
bounds on the local energies W{x; R), Eq.(3.10). For 
instance the speed \q\ of a particle q,q in x is bounded 
by < vfW{x, q, R), R> 0, and its distance p from the 
walls hy p> ^W{x,q,R)-i, R>0. 

In the geometries described in Fig.l we shall prelim- 
inarily discuss bounds at time 0, e.g Eq.(5.3) below, 
and then we shall use energy conservation to extend the 
bounds to positive time. 

Superstability of the potential implies that the num- 
ber N of points in a region A, a cube or a ball, can be 
bounded in terms of the potential energy U in the same 
region and of cpo = Lp{0),lp = ip{^) > (defined after 
Eq.(2.1)). The bound derivation is recalled in Appendix 
A for clarity, and yields the inequality 



Na < C- 



w 



c 



ief 



(5.1) 



,def 



Calling £ = £i/d{x), Eq.(3.11), consider the sequence 
of balls A„ = B{0, 2"), of radii L„ = T^r^, n > nQ, see 
paragraph following Eq.(2.1). Given a configuration x let 
N{x; ^, R) the number of particles in the ball of radius 
Rr^ centered at ^ and 



(1) Vn = { max velocity in A„/f2o )/wi 

(2) pn = min distance to d{Qj n A„) of qi ^ (5.2) 

(3) Ml = max N{x;qi,l) 

(jiGA„/f2o 

Such quantities can be bounded in terms of the maximum 
of W{x, ^, R) over ^ € A„ and, by choosing R = n^/*^, we 
obtain, see definitions Eq.(3.8), 



y„ = {n£p, p„ = {n£)-^,J\fn < C{n£)^ (5.3) 

under the assumption that the wall potential has range 
and is given by Eq.(2.2); the last inequality is a con- 
sequence of the definition of W and of the above men- 
tioned superstability, Eq.(5.1). By controlling the growth 
in time of the energies W we shall extend the validity of 
Eq.(5.3) to positive times. 

Constants convention: From now on we shall encounter 
various constants that are all computable in terms of the 
data of the problem (geometry, mass, potentials, den- 
sities, temperatures and the (arbitrarily) prefixed time 
6): Eq.(5.1) gives a simple example of a computation of 
a constant. To avoid proliferation of labels all constants 
will be positive and denoted C,C' ,C" , . . . , B, B' , . . . or 



8 



VI: Infinite volume. Prictionless dynamics 



c,c' ,c" , . . . ,b,b' ,b" , . . .: they have to be regarded as func- 
tions of the order of appearance, non decreasing the 
ones denoted by capital letters and non increasing the 
ones with lower case letters; furthermore the constants 
C, . . . , c, . . . may also depend on the parameters that we 
shall name £ or E and will be again monotonic non de- 
creasing or non increasing, respectively, as functions of 
the order of appearance and oi £ or E. 

In this and in the next section we shall study only the 

open regularization x^'^'^\t) s\ '°''x, t < Q, Q be- 
ing an arbitrarily fixed positive time. Therefore ^("'^^ (r) 
will be a finite configuration of particles which are not 

necessarily in A„ (unless t = 0). 

Call Pn{t),yn{'t)TJ^n{t) the quantities in Eq.(5.2) eval- 
uated for s("'°)(t) and let p„(f), V n{t) M n(f) be the cor- 
responding quantities defined by taking the maxima and 
minima in the full A„ (i.e. not excluding points in fio). 

For times Q < t < Q consider An-regularized motions 
(see Sec. Ill) evolving from x with n fixed (see comment 
(6), p. 3). Define 



(2) the inequality holds Vd's. The bounds in Eq.(5.5), 
(5.6) will be repeatedly used. 

By the definition Eq.(3.10) of W it follows that 

Vn{s) < C Rnis)''/^ (5.7) 

Therefore for = 1,2, going back to Eq.(5.4) and solv- 
ing it, i?„(i) is bounded proportionally to R„{0) = n^l"^: 
hence i?„(t) < Cn'l'^ for < i < 6, ci = 1, 2 and 
W < Cn hence V„(t) < Cn^^"^ and the bounds in the 
first line of Eq.(5.8), in theorem 3 below, are immediate 
consequences. 

Theorem 3: It is p„(t) > p„(t),K(i) < Vn{t),Afn{t) < 
77n{t) and up to time t < Q the following inequalities 
hold for d= 1,2 or d = 3 respectively: 

Vn{t)<Cni, 7J'n{t)<Cni, pjt)>cn-i 

— 1 3 3 (5-8) 

Vn{t)<Cn^, 77n{t)<Cn^, p^{t)>cn~^ 



R^{t)'^'n^+ fvnisf-^, (5.4) 
Jo 

where i?„(0) = .gi/d(2") = n^''^ for d = 1, 2 but i?„(0) = 
n^/^ for d = 3, and v\ Vn{s) is the maximum speed that 
a moving particle inside a thermostat can acquire in the 
time interval [0, .s] under the A„-regularized evolution: 
formally viVn{s) = maxq.^no,s'<s |ft(s')l- 

By the choice of ^ it is Rn{Q) > n^/'^ for d = 
1,2,3: hence it will be possible to claim if d = 1,2 
that H^(x(0);C,i?«(0)) < n£{x{0)) < oo with ^q- 
probability 1, see Eq. (3. 10), (3. 11) and Appendix B. If 
d = 3 the somewhat weaker bound W{x{0);^,Rn{0)) < 
n3/2£:(x(0)) < 00 will hold. 

The dimensionless quantity Rn{t) will also provide a 
convenient upper bound to the maximal distance a mov- 
ing particle inside any thermostat can travel during time 
t, in units of r^, following the A„ regularized motion. 

Then the following a priori energy bound, proved in 
Appendix C, holds: 

Theorem 2: The An-regularized frictionless dynamics 

def 

satisfies, for t < Q and if W{x, R) = supj W{x; ^, R): 

W{x^"^'>\t),R^it)) < CRnitf (5.5) 
for n large enough and C > (depending only on £). 

Remarks: (1) hence we obtain also a bound on the force 
that can be exercised on a particle by the others including 
the forces due to the particles in fio (as the latter force 
is bounded proportionally to A^O) hence by a constant): 

Y.\F,^\<CRn{tY'^ (5.6) 

3 



The d = 3 case is more delicate because the inequality 
obtained from Eq.(5.4) using Eq.(5.7) gives a blow-up in 
a finite time. 

Following [8] we shall prove in Appendices D,E (theo- 
rem 7) that Eq.(5.7) can be improved to Vn{t) < CRn(t), 
hence Vn{t), Rn{t) < Cn^/"^, which implies the second of 
Eq.(5.8) for Vn{t),Nn{t), Pn{t)- The distinction between 
p„(t), y„(i),A/'„(t) and 14(t),7V„(i), p„(i) is only really 
necessary in the d~3 case. 

For d = 1,2, the bounds on Vn{t),Nn{t), Pn{t) in the 
first line Eq.(5.8) have been just discussed as a corol- 
lary of theorem 2. For d = 3 the bounds in Eq.(5.8) 
will be proved for Vn{t),Mn{t), pn{t) first and then they 
will be shown to imply the bounds in Eq.(5.8) for the 
/3„(t), y„(i),!A7'„(f): see Appendix E, paragraphs after 
Eq.(E.3) 

VI. INFINITE VOLUME. FRICTIONLESS 
DYNAMICS 

It will now be checked that the n ^ oo limit motion 
exists in the sense of the local dynamics assumption, i.e. 

existence of x'^^\t)'^= lim„^oo x^^'^\t) and a suitable 
form of its uniqueness. 

The equation of motion, for a particle in the j'-th con- 
tainer (say), can be written as 

g("^°)(t) =ft(0)+fft(0)+ [\t-T)fi{x^^^°\T))dT (6.1) 

Jo 

where the label j on the coordinates (indicating the con- 
tainer) is omitted and fi is the force acting on the se- 
lected particle divided by its mass (for j = it includes 
the stirring force). 



VII: Entropy bound. Thermostatted dynamics 



9 



Existence of the dynamics in the frictionless, open 
boundary case will be discussed proving 

Theorem 4: If d < 3 and x € 'Hi/d the thermodynamic 

limit evolution x^*^\t)i = limn^ocX^"'^\t)i exists. 

Proof: (adapted from the proof of theorem 2.1 in [8, p. 32], 
which applies essentially unaltered). Let 



Si{t,n) 



def ,_(n,0) 



^(n+1,0) 



(6.2) 



then Eq.(6.1) yields 



5i{t,n) < f -dT {F^S,iT,n) 

J' (6.3) 

3 

where F^^ = boimds the maximum gradient of 

the walls plus the stirring forces (see Eq.(5.8)) for d = 
1, 2 and for d = 3 we can take = ni^; F' = 

maxq\d^ ip{q)\; and the sum is over the number A7"„ of 

the particles 9j(t) that can interact with ^j(r) at time 
T. The latter, by Eq.(5.8), is AT^ < Cn^/^ (rf = 1,2) or 
< Cn3/4 {d = 3) for both s("'0)(t) and x^''+'^'°\t). Let 



V=^d + -), 2'=^1i/2'=+r„ (6.4) 
2 a 

where r„r,^ is the maximum distance a particle can travel 
in time < 6, bounded by Eq.(5.8) by C r^ n^/^ (for d = 2 
the T] could be taken = (1 + -)). Then 



Uk{t,n) 



< C 



Jo ''"'P 



n) ds 



(6.5) 



[C is a function of £ as agreed in Sec.V). Eq.(6.5) can 
be iterated t times; choosing £ so that 2^ + C In^^'^ < 2", 
i.e. £ = 2Cn^/2 which is £ > c2"/^ 5k<n for n large. 



By Eq.(5.8) Un{t,n) is < Cn^/'^, so that for n > k, 



(6.6) 



for suitable C ^C^c > (n-independent functions of £). 
Hence the evolutions locally (i.e. inside the ball A;;) be- 
come closer and closer as the regularization is removed 
(i.e. as n — > oo) and very fast so. 
If qi{Q) e Afc, for n > it is 



showing the existence of the dynamics in the thermody- 
namic limit because also the inequality, for n > fc. 



< C2-2 



(6.8) 



follows from Eq.(6.6) and from q'i^''^\t) — q'i^^^''^\t) = 
lo (/i(9^"'°Hr)) - fi{q^''+'^'\r)))dT. Or, forn > fc and 
|(Zi(0)| <r^2^ 



m <Ce 



-c2 



,n/2 



(6.9) 



def 

calling \xi - x'^] = |ft - qi\/vi + |ft - qi\/r^. 

Hence existence of the thermodynamic limit dynamics 
in the frictionless case with open boimdary conditions is 
complete and it yields concrete bounds as well. Unique- 
ness follows from Eq.(6.3): we skip details, [13]. Hence 
we have obtained: 

Theorem 5: There are C{£),c{£)~^ , increasing func- 
tions of £, such that the frictionless evolution satisfies 
the local dynamics property and if qi{0) G Afe 



|t^"'°^(t)| <«iC(5)/e5, 

distance(# '"^t), a(U,f2j n A)) > c{£)k-^r^ , , ^ , 

(d.IU) 

M(t,n) < C(£:)A;3/4 
\xr\t)-xr\t)\<Ci£Ke-<^^^^''' 

yn>k. The x^°^(t) is the unique solution of the fric- 
tionless equations satisfying the first three of Eq.(6.10). 

For d= 1,2 in the above proof the bounds in the first 
line of Eq.(5.8) could be used instead of the weaker ones 
in the second line (which holds also for d = 3): in this 
way the exponents ^ and | can be improved to ^ and 
i respectively. 

It would also be possible to show the stronger result 
that s(°)(t) e -Hi/d for d = 1,2 and x^°\t) € -^3/2 for 
d = 3: but, for the proof of theorem 1, theorems 2-5 
are sufficient, hence the proof of the stronger property is 
relegated to theorem 11 in Appendix F. 

The corresponding proof for the thermostatted evo- 
lution will be somewhat more delicate: and it will be 
weaker as it will not hold under the only assumption that 
£{x) < 00 but it will be necessary to restrict further the 
initial data to a subset of the phase space (which however 
will still have /Uo-probability 1). 



VII. ENTROPY BOUND. THERMOSTATTED 
DYNAMICS 



00 
n=k 



it) 



-t^'^W) (6.7) 



In this section the regularized motion, thermostatted 
and with elastic boundary conditions, will be compared 
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to the thermodynamic hmit of the frictionless dynamics 
with open boundary conditions whose existence and main 
properties have been established in Sec. VI. 

A direct comparison of the two evolutions along the 
lines of Sec. VI does not seem possible because the quan- 
tity aj is not small enough compared to the "Lyapunov 
exponent" rf in Eq.(6.5): the extra aj introduces a small 
inhomogeneous term in attempting a derivation of the 
analogue of Eq.(6.5). The extra term is of order 2""'' 
which is amplified at rate Cri^, 77 > 1 by Eq.(6.4), to a 
quantity 0(2~"''e'-^®"'') diverging with n. 

The strategy will be to introduce a stopping time Tn{x) 
in the A„ regularized dynamics defined by the time in 
which either the total kinetic energy in any thermostat 
or a local energy reaches a conveniently fixed threshold. 
The threshold will be defined so that within T„(.t) the Uj 
are small and the frictionless and thermostatted dynam- 
ics are very close at least well inside A„. We shall then 
conclude that the number of particles M and their max- 
imal speed Vvi in a region T) of width r^p adjacent to rig 
are, in the thermostatted dynamics, bounded in terms of 
the respective values in the frictionless motion with the 
same initial conditions (finite as shown in Theorem 3 for 
S{x) < E). See Eq.(7.3), (7.4) below. 

Thus before the stopping time Tn{x) the entropy pro- 
duction in the thermostatted motion is simply controlled 
(by Eq.(3.4) it is < CMV), showing that the measure 
/Uo is "quasi invariant". This will imply that with large 
Ho probability the thermostatted motion cannot have a 
too short stopping time and, actually, that the stopping 
time cannot be reached before &. The precise statements 
can be found in the proof of theorem 6. 

Restricting attention to the set Xe C Hi/ 4 of ini- 
tial data Xe {x\ £{x) < E;^ the constants C,C',. . ., 
c,c', . . . will be functions of E as stated in Sec.V. 

Let piio iO be the distance of ^ ^ fio from the boundary 
d^lo and let 



A. =^{^ : PnoiO < r^}, A« : Pn^iO < (7-1) 

By the result in Sec. VI there are M and V (which depend 
on E) so that for all x € Xe, with the notations Eq.(3.9), 
for n large enough: 



maxTVA. ^(^^"'"^a;) <M 
t<e 

max VA..(^|"'°^a;) <V-l 



(7.2) 



Fix, once and for all, k > smaller than the mini- 
mum of the kinetic energy densities of the initial x in the 
v thermostats (x- independent with /io-probability 1, by 
the "no phase transition" assumption, i.e. by ergodicity, 
see comment (c) p. 4). 

Let C| the cube with side r^, centered at a point ^ in 
the lattice Ti^Z'', and using the definitions in Eq.(3.9) 



II II def max(iVc (a;),£Q(x)) 

\m\n = max — -—— , (7.3) 



where ec^ix) = y^ec^Jx). Fixed 7 once and for all, ar- 
bitrarily with i < 7 < 1, define the stopping times 



Tnix) =^ max {t : t < Q : W t < t 

Kj^niSr ' ^ x) 



(7.4) 



>«2"^ \\si"''^x\\n<{logny}. 



Theorem 6: If d < 3, 3 C,C' ,c> depending only on 

E such that for all n large enough: 

(1) For all X G Xe, t < Tn{x) and all gi(0) S A(iog„)T. 



(7.5) 



furthermore with the notations Eq.(3.9), for n large 
enough and for all t < Tn{x): 



7Va.(5["'''x) < M, VK^Sr^'x) < V (7.6) 



def 

(2) the fio -probability of the set B = {x\x € Xe and 

Tn{x) < 9} Z5 



Mo 



(B) < (7e"'=('°en)^'^-FC"My_ 



(7.7) 



Assuming theorem 6 it is immediate to see that it im- 
plies the assumptions of theorem 1 (locality of the dy- 
namics) because 7 > 1/2 in Eq.(7.7) will imply (by 
Borel-Cantelli's lemma) that, with /xo~probability 1, 
eventually for large n, Tn{x) > 6 and therefore, by 
Eq.(7.5), the thermodynamic limits of a;(">'')(t), 
will coincide for t < 6. Hence theorem 6 implies theorem 
1. 

Remarks: (1) The proof of theorem 6 contains a standard 
part, namely proving Eq.(7.5) and its corollary Eq.(7.6) 
for t < Tn{x): it follows from the existing literature, see 
[5, 8] and the proof in Section VI, and for completeness 
it is discussed, providing the details, in appendix K. 

(2) The proof of Eq.(7.7), i.e. the proof that the Eq.(7.6) 
is not an empty statement, is the technically original con- 
tribution of this paper: It is based on an "entropy es- 
timate" and the strategy of its use is discussed below, 
deferring to the Appendix J the actual computations to 
obtain the estimate. 

(3) The proof below refers only to the case o = 1: the 
same argument applies for the (easier case of) frictionless 
motions with elastic boundary conditions, a = 0. 
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Proof: The stopping time T„(x) is the time when the tra- 
jectory s[^'^^x, X G Xei crosses the (piecewise smooth) 
surface S' of points y where cither the kinetic energy in 
some ilj n A„ has the value k2"'' or crosses 
from below the value (logn)'''. Denote by S the subset 
of S' of all points 5'^"'(t\a;, x e Xe, T„(a;) < 6. Thus 

E = U^<eS'r"'^^S^ where 

i:/^ {x&XE\Tn{x)=T) (7.8) 

and B is the disjoint union of over r < 6. 

The surfaces S, S' and St, r < ©, are symbolically 
described in Fig.2. 




Fig.2: The horizontal "line" represents E', E is the subset 
of E' under the "curve" S^^g^^-^x; 0{x) < Q is such that 

^^-(kI)^ G Xe while S'l^'-^x ^ Xe for all t € (6'(x),e]; the 
vertical line represents the trajectory of the point S^^^^^^x, the 
incomplete ("dashed") lines the "levels" Et-jE^/, their miss- 
ing parts are made of points not in Xe hut with an "ancestor" 
in Xe- 

This is the setup described in Appendix H, with E a 
"base" and 9{x) = maxT-<e {t\s'^^^^\ G St} a "ceil- 
ing function", so that (see Appendix H and, for early 
applications of this key estimate, see [5, 6, 14]): 



Mo(S) < / yo,i:{dy) / dt w{y) e"(^'*) (7.9) 
Jn Jo 

where (1) i/o,s' denotes the surface area measure (mea- 
suring the volume with /io) on S', 

(2) w = \vx-nx\ with Vx the i^"'^) evaluated at x (by the 
equations of motion) and is the normal to S at a;, 

(3) d{y,t)'^= /°Ja(y("'i)(t'))|rft' bounds the phase 
space contraction, i.e. the entropy production. 

Effectively this means that the distribution jjLQ can 
be treated as an invariant one for the purpose of es- 
timating probabilities in B via Eq.(7.9), because phase 
space contraction a = a{x) is given by Eq.(3.3) and Qo, 
see comment (2) following Eq.(3.3), can be estimated 
in the same way as Q^, by C MV (i.e by a bound on 
the speed times the number of particles in A» times the 
number of particles in fio)- It follows that the integral 
j-T„(a:) |^^2;("'")(t))| is also uniformly bounded (in B) 
by CMV. Therefore a volume element in B contracts at 
most by g"*^^^ on the trajectory of /xo-almost all points 
a; G B, up to the stopping time Tn{x), and it follows 



Mo(S) < e^'^''^^, \w{y)\vQ^^,{dy) (7.10) 

which we call, naturally, entropy hound. 

The derivation of the Eq.(7.7) becomes, by Eq.(7.10), 
an "equilibrium estimate" , as it does not involve times 
t > and is a standard consequence of the superstability 
property of the potential ^p: the detailed calculations are 
in Appendix J. 

Remark: The analysis in this section is valid for all di- 
mensions except for the use of the properties of x^'^^(t), 
theorem 6, which have been derived only for d<2>. 

VIII. CONCLUSIONS 

Equivalence between different thermostats is widely 
studied in the literature and the basic ideas, extended 
here, were laid down in [15]. A clear understanding of 
the problem was already set up in comparing isokinetic, 
isoenergetic and Nose- Hoover bulk thermostats in [15], 
where a history of the earlier results is presented as well, 
see also [1, 16, 17]. 

There are, since a long time, studies of systems with 
free thermostats, starting with [2]. Such thermostats are 
somewhat pathological and may not always lead to the 
stationary states that would be expected: as exemplified 
in the case of simple spin chains, [18, 19]. More recently 
similar or identical thermostat models built with free sys- 
tems have been considered starting with [20] . 

Isokinetic thermostats should be treated in a very sim- 
ilar way, [1] : the extra difficulty is that the entropy pro- 
duction in a finite time interval receives a contribution 
also from the time derivative of the total energy of the 
reservoirs, [1], and further work seems needed. 

More general cases, like; Lcnnard-.Joncis interparticle 
potentials are difiicult, see [21], and new ideas may be 
needed. 

Finally here the interaction potential has been assumed 
smooth: singularities like hard core could be also consid- 
ered at a heuristic level. It seems that in presence of 
hard cores plus smooth repulsive potentials all estimates 
of Sec.V are still valid but the existence of the limiting 
motion as A ^ oo remains a difficult point because of 
the discontinuities in the velocities due to collisions. 

It should be noted that the key bounds Eq.(5.8) hold 
for the open boundary conditions motions x'^^'^\t) (as 
needed for our purposes) and have not been proved for 
x^^'°'\t). A careful analysis of our argument shows that 
the bounds hold only for particles initially in A^, k < n, 
in the cases of reffecting boundary conditions, but not 
for k = n: the results of Sec. VII show that the maxi- 
mal speed Vn{t) is boimded by C n^/^(log n)''' in the A„- 
regularized thermostatted or frictionless dynamics rather 
than by Cn^/"^, which is the bound obtained for the fric- 
tionless dynamics with open boundary conditions. 
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The main problem left open is what can be said about 
the limit t — > oo, i.e the study of the stationary state 
reached at infinite time. A conjecture has been proposed, 
[22], that can be interpreted as saying that if d = 1,2 
the limit will be "trivial" : i.e it will be an equilibrium 
Gibbs distribution at some intermediate temperature if 
# = 0. But, again interpreting the conjecture in [22], 
for d = 3 the stationary distribution (s) will be nontrivial 
and asymptotically they will be Gibbs distributions at 
the initial temperatures and densities of the thermostats. 



Appendix A: Superstability. Sets of full measures 

In our simple case {ip > 0, decreasing and < oo) 
the superstability bound is reduced to Schwartz' inequal- 
ity. In fact, by the dimcnsionless energy definition in 
Eq.(3.10), W>{^+N)>^ with the sum run- 

ning over labels p of disjoint boxes of diameter ^ cover- 
ing a cube or a ball A and containing Np > particles (in 
particular: N = Y^p^p), hence over i < \A\{2Vd/r^f 

terms. By the Schwartz' inequality y/Tp/2ipoN < VW £, 
hence Eq.(5.1). 

There are cq and Rq and a strictly positive, non de- 
creasing function 7(c), c > cq, so that Vc > cq,\IR > Rq, 



Appendix C: Theorems 2 (d < 3) and 3 (d < 2) 



Consider the ball i?„(f, s)) around ^ e M'' of radius 



R„{t,s)=^Rn{t)+ r^^^^ds>l. (C.l) 



The ball radius shrinks as ,s increases between and t at 
speed vi 14, (s): therefore no particle can enter it. 

The quantity Rn{t, s) can be used to obtain a bound 
on the size of TV(a;("'°^ {t); ^, -Rn(t)) in terms of the initial 
data x{0) = X and of 



W{x,R)'^= supW{x;^,R). 



(C.2) 



Remark that W{x^"'^\t), R) is finite if the just defined 
meaning of x^"''°\t) is kept in mind (so that x^"'°^(t) 
contains finitely many particles). 

Lot x^iQjR) be a smooth function of g — ^ that has 
value 1 in the ball B{^, R) and decreases radially to reach 
outside the ball B{^,2R) with gradient bounded by 
2{r^R)-'^. Let also 



l,J{W{x,Q,R) > cR'^}) < e-^(^)^'^ (A.l) ^^^5^'^)'- ^ E ^d<l,R) 



(which more generally holds for superstable but not 
necesssarily positive cp, [9]). 

If 5 : Z'* — ^ M+, g{i) > 1, c > co, the probability 

^'o{^ieZ',reZ,r>,ii) {W{x;i,r) < cr^}) (A.2) 

is > 1 - E,eZ^^eZ,^>g(i) ^^^^ ^^^^S 

bounded proportionally to the sum X^^g^'' e"'''^'''!^^''! 

which converges if g{i) > c'(log_|_ |i|)"^/'', with c' large 
enough. 



, mq^ , , , 1 
(— +V'(g)+2 



. def 



E 



fiQ - 1') + (Po), 



W^(x;i?) = supVF(a;;^,i?). 

(C.3) 

Denoting B an estimate of how many balls of radius 1 
are needed in R'' to cover a ball of radius 3 (a multiple 
of the radius large enough for later use in Eq.(C.7)) so 
that every pair of points at distance < 1 is inside at least 
one of the covering balls, it follows that W{x;S,,2R) < 
BW{x,R), see Eq.(C.2), so that: 



Appendix B: Choice of Rn{t) 

The proof of the inequalities Eq.(5.5) yields \/t < Q 
that W{sf'°\,R) < cW{x,R+f^ Vn{T)dT/Q) provided 

R+ j* VndT/e 



R is such that 



< 2, which is implied by i? > 



i?o + Jo Vn{s)ds/Q, Rq > 0. The maximal speed viVn{t) 

at time t is bounded by V„{t) < \/w(s["'°^x,R). The 

choice Rq = Rn (0) = ni/'* is the weakest that still insures 
that the set of initial data has W{x; 0, R)/R'^ finite with 
/UQ-probability 1, see Appendix A. 



W^(x("'°)(T);^,i?) < VF(s("'°HT);^,i?) 

< W(s("'0)(r);e,2ii), 
W{x\i,R) < BW{x,R) 



(C.4) 



Although W has a direct physical interpretation W turns 

out to be mathematically more convenient and equiva- 
lent, for our purposes, because of Eq.(C.4). 

The proof of theorem 2 is taken, from the version in [8, 
p. 34] of an idea in [5, p. 72]: it is repeated for complete- 
ness because quite a few minor modifications are needed 
here. It could be extended, essentially unaltered to the 
other boundary conditions, see [23, Sec. 6- A], but we shall 
not need it here. 

Consider W(xa„(s); ^, Rn{t, s)), for < s < i < 6: 
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q'mo 

(C.5) 

because the s-derivative of X5(9(s): Rn{t, s)) is < since 
no particle can enter the shrinking ball B{^^, Rn{t, s)) as 
s grows: i.e. XiiQi^)', Rn{t, s)) cannot increase with s. 
The sums are restricted to the q's of xa„(s)- 

In the frictionless case a computation of the deriva- 
tive in Eq.(C.5) leads, with the help of the equations of 
motion and setting xdl{s),Rn{t, «)) = Xi,q,t,s, to 



q{s)dqip{q{s) - q\s)) 



ds 



2(^0 



(C.6) 

where F{x{s)) denotes the force that the particles in J7o 
exercise on the thermostats particles and the dot indi- 
cates a s-derivative; keep in mind that positions and ve- 
locities of particles outside A„ are not considered, in the 
A„-regularized dynamics at open boundary conditions. 

Since the non zero terms have \q{s) — q'{s)\ < r^p, 
the gradient of % is < 2 {r^pRn(t, s))^^ and |(?|, Ig'l < 
viVn{s) = rip\Rn{t, s)\ it follows, since F — max(A'o|i9(p|) 
is an upper bound on the force that particles in can 
exercise on a particle in the thermostats, 



d — Fvi — 

as ipo 

Fr^ l^"(^'")I.BW^(xA„(a);^,2i?„(t,s) + 1) (C.7) 



+ 



ifO Rn{t,s) 



< £2 Fvi \Rn{t,s)\ 

~ ipo V\ Rn{t,s) 



+ l)W{XA^{s);Rn{t, s)) 



where W{x;R) is defined in analogy with Eq.(C.2), and 
the first term in the r.h.s. is obtained by bounding 
the corresponding term in Eq.(C.6) as J2q{\Q\\/X)\/X fol- 
lowed by Schwartz' inequality. 

By Eq.(C^),(C.4), R„{t,s)/R{t,0) < 2, VF(xa„(s),^, 
Rn{t,s)) < W{x\^{s),Rn{t,s)) ws get the inequality 

W(SA„(.s),i?„(t,s)) < CW{xKM.Rn{t,0)) (C.8) 



withCt/exp((^log2+3|i^)eB). 
This concludes the proof of theorem 2. 

To prove theorem 3 for d = 2 the inequalities Eq.(5.8) 
for pn{t),Vn{t),Afn{t) follow from Rn{t) < C'n}''^ and 
from Eq.(5.5),(5.7). It remains to check them for 



Pnit), Vn{t),JJ'n{t), i.e. it remains to consider the parti- 
cles moving inside Oq. 

The inequalities certainly extend to the initial time t = 
and to any t for A7'„(t) because the number of particles 
in J7o is fixed and bounded in terms of £. The variation 
E{t) — E{0) of the energy contained inside the region fig 
equals the work Cext done by the particles outside 
plus the work of the stirring forces in the time t. By 
the bounds Eq.(5.6),(5.7) this implies, for t < Q: 



(Afo|l$||niax|g(r)| +Aroni/2(7y„(r))dr (C.9) 

qeflo 



Therefore if v{t) = ma.Xq^a„ ,t' <t \q{t')\ if follows that 
vit)"^ < C'n + C"v{t) which implies v{t) < viCn^/'^; 
hence il>{q{t) < Cn and the Eq.(5.8) are checked for 
d = 2. 

Theorem 3 will be complete after the proof of the sec- 
ond line of Eq.(5.8) which follows from the speed bounds 
in Appendix E. 



Appendix D: Work bounds. Thermostats in d = 3 



In this and and in the next appendix we study the 3-d 



dynamics S^^'^^x 



,(n,0)^,) ^ {t^'\t),qr'{t)), fne- 



tionless and with open boundary conditions. We follow 
the method of [8]. 

Time intervals will be often subdivided into subinter- 
vals of equal length 6: in general, however, the number 

may be not an integer. We shall intend, without mention, 
that the last of the intervals might be shorter. This is 
used, without mention, only in cases in which the bounds 

considered remain valid also for th e "left over" interval. 
Let e{q, q) '^^ ^j^^{fq^ + tP{q)), abridged as e^it) =^ 

e{qj {t),qj (t)) when q, q are time dependent. Remark that 
the speed of a particle is < vie{q, q) and its energy in the 

2 

field of the walls forces is ^^^e(<7, q)'^. 

A particle qi ^ fig will be called fast in a time interval 

def 

JC [0,t] if min^rsJ > Rn{t)- 

The main property on which the analysis rests is an 

estimate by CRn{tY~^ of the work C'^^ Ij^j 'ijFij 
that the i-th particle, assumed "fast" , performs over the 
neighbors in the "short" time interval J = [s — 6,s], with 

O'^if !j£i?„(t)-7 with s < t and 1 > 7 > 0, [8]. And the 
precise property that we shall need is expressed by the 
lemma: 

Lemma 1: Consider the dynamics s["'^^ in the geometry 
of Fig.l and in a time interval J C [0,©] of size 6 = 
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^Rn{t)-^. Then the work C = Jj J2j QjFij that the i- 
th particle, assumed "fast" and in flj, j > 0, performs 
over the neighbors in the time interval J is bounded by 
m < CRnitf-^ if -re (i, |] and J c [0,t]. 

This will be proved considering the special geometry in 
Fig.l, following [8]. 

All things considered it will turn out that if d = 3 the 
parameter J can be arbitrarily fixed in the interval |]. 

Remark: The restriction on the geometry is severe: it is 
used only to prove lemma 5 below. It is very likely, how- 
ever, that lemma 5 holds for all thermostats considered 
in Sec. I, and that our equivalence holds in the general 
cases as well. 

Notice that all particles Qj interacting with a given 
particle at a given time r are contained inside a ball of 
radius 2r^Rn{T). 

Identifying particles with their labels and imagining 
the following maxima to be taken over t € J, divide the 
particles {qjAj) into three groups, [8], as "slow", "inter- 
mediate" , "quick" : 

(a) Ao = slow — {j| max |£j(t)| < Rn{ty}, 

(b) Ak = intermediate = {j\2^~^Rn{ty < max|ej(T)| < 
2'^Rn{t)'^}, for = 1, . . . , fcm — 1, with km being 
determined by the condition 2^'^Rn{t)'^ = 

(c) A= quick = {j\ max \ej{T)\ > ^Rn{t)}, 

Correspondingly the work C = Cq + Ck + JC: the 
three contributions will be estimated separately. 

This corresponds to proposition 4.3 in [8]. The time 
interval J = [s — 9, s] (s > 9) is assumed, for simplicity, to 
be of length 9 = Rn{t)~^r^/vi. The estimates, however, 
hold also for shorter intervals. 

C-1 Work Co of fast particles on slow particles 

Divide the time interval J into H consecutive intervals 
A/j = [sh,sh+i\ of size equal to 2viRn(i) ^^"^ ^ ~ 
2Rl^^ be their number. The work Cq is bounded by 

l-Col <\\d^\\viRn{ty j N{T)dT 

H ' (D.l) 
<Ci?„(ir V|A^|maxJV(T) 

^-^ reAh 
h—1 

where: N{t) is the number of slow particles in the ball 

of radius r^p and center g^"'"'' (r). The time intervals of 
size 9 = Rn{t)~^r^/vi are so short that a slow particle 
cannot travel a distance greater than the range away 
from qj{s — 9). 

A slow particle contributes to N{t) if |5j"'°^(t) — 
#'°^(r)| < r^, hence \qf^'\s -6)- #'°)(t)| < 2r^ 
because — 6) — 9j"'°^(T)| < r,^. 



If Nh = max^eAh N(t) then: 

|£o| < CRn{t)-^^-^y\Y. ^hY'^ 
fe 

having used Ef=i Nh < sfH^f^J^ and H = 2R\-'i. 

Notice that Nl is bounded by W{x{s-e)). qf'^\T).2) 
where r is the time when the maximum defining Nh is 
reached. 

Hence, if Th denotes the tube spanned by the ball of ra- 
dius 2r<^ and center at g (i) for t € A^, is bounded 
by W{rh) provided W{T) is defined as in Eq.(3.10) with 
B{^,R) replaced by T. 

Lemma 2: Each Th can be intersected by at most C 

other Th' ■ 

Remark: Thus Tfi is a set contained in a ball of ra- 
dius ry,Rn{t) (which is, by definition, the maximum 
distance any particle can travel) and each point is in 
at most C different Th's. This implies ^f^W{Th) < 

CW{T) < CERnitf, if r =^ Ufe Tft. Concluding: 
{Y.h^hf'^ < CW{T)^'^ so that by Eq.(D.2) it is 

I'Col < CifaRn ^ R^/"^, d = 3, because T is contained 
in the ball of radius Rn{t) r^ centered at g-"'°''(s — 9), or 

\Co\ < C'foRnitfRn'^'''^^ < C^oRnitf-^ (D.3) 

by Eq.(5.5), provided 1 — ^ > 7- This is possible for 

rf = 3 if 7 < |. 

Proof: To check lemma 2 let r — >■ Ao(t) be the path that 
the fast particle would follow under its own inertia and 
the walls forces, which act in a small strip V of width 
< r^ near the containers walls, in a time interval of 
size 9 = Rnit)~'^r^/vi. 

Given a point ^ e the particle following Aq (t) might 
spend up to two intervals of time Ji , J2 C J inside the 
ball of radius r,^ centered on ^ (possibly one before the 
collision with the wall and one after bouncing) . However 
the actual path r ^ A(r) differs from Aq. 

Suppose that initially the fast particle is outside 2?: 
until outside V it proceeds, in time < 9, over a distance 

e{T)>Msh)\{r-sn)- r{T-s)^-^ds 

JsH ^ (D.4) 

>|t - Sh\ (t^min - CvrRn{t)-'<+'^/^), 

by Eq.(5.6), undergoing a deflection by an angle which 
is at most CRn{t)~^~^'''^'^ / Rn{t) and will have velocity 
> ^Rn{t), for n large. 

At the entrance into the region 2? at a time ti the 
difference between X{ti) and Ao(ti) will be < \ti — 
Sh\C viRn{t)~^~^'^^'^ and there are two possibilities for 
n large. The component of the velocity along the normal 
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to the wall is > -^f i or its tangential component 

• \ 1 Rn{t) 

In the first case the time spent by Ao inside T> will 
be < ■^^^p^ Rn{t)^^ and then at a time t2 the particle 
will be again out of V with a position |A(t2) — Ao(t2)| < 
(^2 — Sh) C vi i?„(t)~^+''/^ and speed normal to the wall 
> cRn{t). The motion will then proceed a distance 
bounded as in Eq.(D.4) (with a larger C). 

In the second case; the tangential velocity is large and 
the particle will move away from the entrance point into 
X> keeping its velocity component parallel to the wall 
within < CB.n{t)-^'+'^/'^/Rn{t) dashing away from any 
fixed ball of radius r^, without coming close to any point 
within 2rip more than twice. 

Next we can consider the case of a fast particle initially 
in I?: if the velocity component normal to the wall is 
— 2v^^^^"^*^ particle will get out of D in a short 
time < C^Rn{t)~^-, if the tangential component is large 
the particle will move away from initial position in 2?, 
keeping its velocity component parallel to the wall and 
undergoing a deflection < Ci?„ (t) V^n(t)- 

Remark: Lemma 2 is likely to hold under the sole as- 
sumption that the fast particle path Aq moves away from 
the origin at radial speed bounded below proportionally 
to the initial value of . Hence it holds probably in full 
generality for the thermostats considered in Sec. I. 

C-2 Work jCk of fast particles on 
intermediate speed particles 

If N{t) is the number of intermediate speed particles 
within an interaction radius of ^^(t), Ch is bounded by 



\jOk\ < \\dip\\vi2''Rn{ty 



N{t) dT 



(D.5) 



analogously to Eq.(D.l). 

The relative speed of the particles qi and qj , by the def- 
inition of Ak, can be bounded below by |ei(r) — Sj{T)\ > 
Rn{t) — 2'''^~^Rn{t)^ > ^Rn{t) for n large as long as the 
fast particle is outside V. 

The contribution of the j-th particle to Ck lasts, there- 
fore, a time boimded by CRn{t)^^ except for the time, 
also of the order of Rn{t)~^ , spent by the fast parti- 
cle within the wall with tangential speed smaller than a 
fraction of say 3i?„(t)/4. If at tj G J the ej(r) is 

maximum, then 



\Ck\<\Ak\CRnit)-^2''Rr^{ty, and 



(D.6) 



< 



^ (e^is-ef + cj^^ ^dr\qj{T)\N^ir)) 



where Nj (r) is the number of particles of (r) which 



interact with Q'j(r). We multiply both both sides of (D.6) 
by 2-^, bound |gj(T)|2-'= < viRn{t)'^ and 



J2 E ^i(^) ^ '^^"W". d = 3 (D.7) 

k jeAk 

which can be proved as follows: let B',B be the balls 
with some center qi{s — 6) and radius 2Rn{t) and 3ii„(t), 
respectively. Then 



EE^.(-)^ E E 1 



(D.8) 



< 



where the Cj's are elements of a covering of B by disjoint 
cubes of diameter r,^/2 and Nd are the numbers of parti- 
cles in the Cj's; and the ^* denotes sum over Cj, Cj C B 
with dist{Ci, Cj) < r^. The inequality continues as 



c.cB (D-9) 
< W{x{T),qi{s - e),Rn{t)) < CRnitf 

(by Eq.(5.5)) which proves Eq.(D.7). Since J2k 
Y.j£j{T? < CRn{tf, by Eq.(5.5), we get from Eq. 
(D.6)), using Eq.(D.7): 

Ei-^fe|2'=-'^n^<E2-'E-i(«-^)' 

k k j(^Ak (D.IO) 

+ C'Rn(tyRr,{tfRr,{t)-'^ < C" Rn{tf . 

Therefore J2k IA|2'= < C{Rl~^^) < CRn{t)^-^0/Q and 
J2k ^k is bounded, via Eq(D.6), proportionally to 



Rn{tyRn{t)-^Rn{tf-'-^ = CRn{tf"' (D.ll) 



C-3 Work £ of a fast particle on quick particles 

The work (of the inner forces) on the quick parti- 
cles can be bounded by first remarking that their £j(r)^ 
at a time r e J is bounded below by ^i?„(f)^ minus 

CRn{t)i~'^Sj (upper bound on the work done by the 
other particles in time ^i?„(t)~^). Therefore it follows 
that ej{T) > cRn{t) for n large enough, since 7 > |. 

So that \A\{cR^{t)Y < |g,(T)P < CvlR^itf, by 
Eq(5.5). _ 

It follows that 1^1 < C£Rn(t). Therefore by the en- 
ergy estimate (see Eq.(5.5)) an application of Schwartz' 
inequality bounds \C\ above by 
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< C'\A\^ / dr( V ej(r)2)i/2 < C Rn{tf^ Rn{t)i 9 

Js-9 

(D.12) 

by Eq.(5.5), i.e. \C\ < CR„{tf-'i. 

Appendix E: Speed bounds in d = 3 



so that, by the work estimates and by H6 < 6, it follows 

G^R^itf < ClRnitf + C" Rnitf (E.2) 

for n (hence Rn{t)) large enough if C" is the largest of all 
constants met so far. Since G > 2Ci, Eq.(E.2) implies 
< |C" and the proof of theorem 7 is complete. 
Going back to the definition of Rn{t) in Eq.(5.4) it 
follows that 3 C > 0, depending only on £, such that for 

t < e, 



We can now infer an estimate for maximal speed, 

viVn{t)i 'iiid maximal distance, i?„(t)r^, that a thermo- 
stat particle can travel in the A„-rcgularizcd motion. 

Theorem 7: Given x G Hi/d, d — 3, there is a constant 
G depending only on £ = £{x) and Q such that for any 
particle of x initially in A„ it is s{t) < G Rn{t) for all 
t < e. Therefore, by Eq.(5.4), Vn{t) < CR„{t) and 
Rn{t) < CRn{0) = Gn'l'^ fort<e. 

Proof. Suppose that there arc f* < t < S and qi{0) E A„ 
so that e{t*) = GRn{t). We claim that the above cannot 
hold if G^ > |C", where G" is a constant defined below, 
in the course of the proof. 

Remark that the claim implies theorem 7 with G = 
(4C")/3)^/^. And the claim is proved as follows. Let Ci 
be such that £,(0) < Cii?„(0) for all x with £{x) < E. 
Then by choosing G" large enough it will be G > 2Ci. 
Since ei{t*) = Gi?„(i) > GiRn{t) > ei(0) there exists 
tx < t* and £,(t) > Gii?„(i) for all t e [ti,t*]. We shall 

— 1/2 

next prove that t* — ti > GRn so that [ti, t*] splits into 
the union of H intervals of length 9 = ^R„{t)~'^ with 

H large. In fact < Gsi{T)^\Fi{x{T))\, hence 

< C ^\F,{x{t))\, and by Eq.(5.6) (for non integer 
H see the comment at the beginning of Appendix D): 

GRn{t) < GiRnit) + G{t* - t,)R,,{tf/^ (E.l) 

which yields t* ~ ti > cRn^ > Rn''^ if 7 > i, d = 3. 

The variation {ei{t*)'^ — ei(ti)^) equals the work of the 
pair interaction forces over the i-th particle. It is there- 
fore the variation of the pair-potentials energy of the par- 
ticle Qi in the configurations x^"'^\ti),x'^"'^\t*) not tak- 
ing into account the A^o particles in flo, plus the work 
that the particles in flo perform on the particle qt, plus 
the work that particle i performs on the neighbors. The 
three terms are bounded by 

(a) \\dip\\ times the number of particles interacting 
with i at the two times, hence by GRn{tY^^ and by 
Eq.(5.6),(5.5) (valid also if = 3, see Sec.V), 

(b) the work performed on the particle qi by the A^o par- 
ticles in Oo is bounded by GNoRnitY^"^ by Eq.(5.5), 

(c) I Eti i::tth-i)0 qj{r)F,j{x{T))dT\ 



W{x^"'°\t),Rn{t))<GRn{tr, 
Rn{t)< Gn'/^, K(t)<Gni/2, 

which means that the maximum speed and the maxi- 
mum distance a particle not in SIq can travel in the A„- 
regularized dynamics grows as a power of n, at most. 

It remains to consider the particles in flo to derive 
the bounds Eq.(5.8) with Vn{t), J\fn{t) , Pn{t) replaced by 
Vn(t),Afnit),'Pn{t). i.c by the quantities in which the 
particles in Q.q are also taken into account. However 
the proof of the corresponding statement for d = 1,2 
(see paragraphs preceding and following Eq.(C.9) in Ap- 
pendix C) applies unaltered as it only depends on bounds 
on Vn{t) proportional to n^^"^. 

It can be remarked that the Eq.(E.3) summarizes re- 
sults that in the cases = 1, 2 were the contents of the- 
orems 3,4 in [23]. However the bound on Rn{t) is sub- 
stantially worse than the corresponding in dimensions 1, 2 
because it only allows us to bound the maximum num- 
ber of particles that can fall inside a ball of radius r,^ by 
Gn'^l^ rather than by Cn^/^ as in Eq.(5.8). 

A key implication of the above results is a proof of 
Theorem 4 in d = 3: 

Theorem 8: If x E 'Hi/d the thermodynamic limit euo- 

lutionx'^^\t)i = \imn-^oo{S^ '^^x)i exists. 

Proof: Remark that in Sec. VI no use has been made of 
entropy bounds. The only difference is that with the new 
bound in Eq.(E.3) r] in Eq.(6.4) is one unit larger. 

By the final remark to Sec. VII the analysis leads to 

Theorem 9: With no -probability 1 in x G 'Hi/d; d = 3; 
the frictionless motions have the local dynamics property 

and, as n —)■ 00, the limits of s["'^^x, S^^'^^x exist and 
coincide, and furthermore the Eq.(2.3) with A = R'' for 
x^°^(t) have a unique solution for data in and with 
values in T-Lz/i- 

This is proved but for the uniqueness statement: as in 
c? = 1, 2 the proof is by iteration of the integral equations 
for two solutions with the initial data in T-Lq and with 
values in "^9^/2: we skip the details, see [13]. 

The analysis in Sec. VII can also be followed word by 
word to check the expected result 
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Theorem 10: The limits lim„_ 



.(n,0) 



X exist and co- 



incide with lim„_ 



-o(".o). 



X. 



The check is left to the reader. 



Appendix F: Frictionless motion is a flow if d = 2 



The following theorem is obtained by a straightforward 
adaptation of theorem 2.2 in [13]. 

Theorem 11: Let E > and d= 1,2. Then, given any 
O there is E' (depending on Q,E) so that for all x such 
that £{x)<E 



Suppose po - pn( > 9 (,{£,/ Tip), hence po > Cn''^, then 
W{xf^;^,po) < C'W{x,po + Pn,) < C"{po+Pn,Y < 
Cpl and M' < Cp'^-^n/^ < C pl~^^^'^'^^\ Hence 

W{xf\tp) <C{pt + pV^'/^''^) < Cpi if c = lid 
for d = 2 and C, = 1/2 if d = 3; likewise one finds C = 1 
if d = 1. 

The values of (n^ — l)*" < po < Cn^ are still to be 
examined. In this case, however, the boimd W{xf^ -.S^, p+ 
Pn/^ < CW{x, p + 2pn^) involves quantities p, pn^ with 
ratios bounded above and below by a constant, hence 
W{xf^;^,p) is bounded by W{x; ^, p + Cnj) < C p"^ . 

Conclusion: there is C > 0, depending only on £ and 
for all p > gd^/r^), i < 6 it is W{x['^^;^,p) < C p'^ if 
C = 1/rf for = 1, 2 and C = 1/2 if d = 3. 



for allt<e 



(F.l) 



so that the evolution x — >■ s[^^x is a flow in Hi/ a- For 
= 3 data in Hi/ 3 evolve into data in "^3/2- Hence the 

same holds for s'f^x for a = 0, 1 hy theorem 7. 

Proof. The proof is based on the bounds Eq.(6.8),(6.9) 
derived in the proofs of the theorems 6,7 (for d = 1,2) 
and 10 (for d = 3). 

Let s["''^^x and consider 



W{x[''\^,p), for p > (log+dCI/r^))^/"^ (F.2) 

with W defined as in Eq.(C.3). Let — 1 be the smallest 
integer such that A„j_i contains the ball of center ^ and 
radius pr^. Then < G 



W{x[°\i,p)<W{xr''>,^,p) 



t 



,0) 



+ \W{x['\^,p)-W{xr^''>,^,p)\ 



(«5,0) 



(F.3) 



The particles which at any time t < <d contribute to 
Eq.(F.3) are in A„^ — 1. Via theorem 6, maximal speed, 
distance to the walls and number of interacting particles 
are bounded by C"2"«/^; hence: 



\W{x\°\^,p)-'w{xr'"\^,p)\<C'e 



(F.4) 



Consider first the case of p large, say p > pn^ ■ Then 

W(x[^^;S,, p) can be estimated by remarking that the ar- 
gument leading to Eq.(5.5) remains unchanged if R(t) = 
p + /o Vn^ {T)dT/e and R{t, s) = R{t) + j' Vn, {T)dT are 
used instead of the corresponding Rn^{t), Rni:{t, s) (as 
long as p > 0). Then 



W{xi°^;^,p + Pr,,) <CW{x,p + 2p^^) 
as in Eq.(5.5) (see also Eq.(C.3)). 



(F.5) 



Appendix G: Free thermostats 



The need for interaction between particles in order 

to have a physically sound thermostat model has been 
stressed in [24, 25] and provides a measure of the impor- 
tance of the problems met above. 

The discussion in Sec. IV is heuristic unless the local 
dynamics assumption is proved. However if the model is 
modified by keeping only the interaction (p between the 
test particles and between test particles and thermostat 
particles, but suppressing interaction between particles 
in the same ^j, j > 0, and, furthermore, replacing the 
wall potentials by an elastic collision rule {i.e. supposing 
elastic collisions with the walls and Uj{^j) = 0, j > 0) 
the analysis can be further pursued and completed. This 
will be referred as the ''free thermostats" model. 

In the frictionless case this is the classical version of the 
frictionless thermostat models that could be completely 
treated in quantum mechanics, [2]. 

Let An be the ball B(0, 2") of radius 2"r^ and n > 
no (see the paragraph following Eq.(2.1)). If N bounds 
the number of particles in the ball B{0,2"°) up to an 
arbitrarily prefixed time O, the first inequality Eq.(4.2) 
and the supposed isoenergetic evolution (which in this 
case is also isokinetic) imply 



\aj\<NoF\ 



N 



< 



NoF 



def , 



2Kj/m ^dkBTj/m 



(G.l) 



It follows that, for ( > 1/d, the speed of the particles 
initially in the shell A„/A„_i with radii 2"r<^, 2"~^^ry, will 
remain within the initial speed up to, at most, a factor 
A^^ = e^^®. The initial speed of the latter particles is 
bounded by, see Eq.(3.11), 



Vn = Vi 



(G.2) 
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Hence, if n(0) is the smallest value of n for which the 

inequality 2"'r^ — y„ A < Dq + does not hold, no 
particle at distance > 2"(^)+^r,^ can interact with the 
test system. 

This means that N < £(; (x) 2 d- and the dynam- 

ics becomes a finitely many particles dynamics 

involving < No + N particles at most. 

From the equations of motion for the Nq + N particles 
we see that their speed will never exceed 



Ve = {V + FNoNe)X 



(G.3) 



if V is the maximum of their initial speeds. In turn this 
means that for n large enough a better bound holds: 



< 



NqNVAF 



s>0 



(G.4) 



with T = minj>o Tj and 5 = minj>o 5j and Wj^„2'^"r^ 
bounds below (for suitable cjj.„) the volume of VLj n A„. 



Hence, for a — 0, 1, it is lim 
and also the dynamics is local in the above sense. This 
completes the analysis of free thermostats and proves: 

Theorem 12: Free isoenergetic and frictionless ther- 
mostats are equivalent in the thermodynanic limit 

Notice that essential use has been made of the prop- 
erty that, in absence of interaction among pairs of ther- 
mostat particles and in presence of perfectly elastic walls, 
isokinetic and isoenergetic dynamics coincide: so the de- 
nominators in Eq.(G.l) are constant. 

It would be possible to consider non rigid walls, mod- 
eled by a soft potential ip diverging near them. We do 
not do so because it will be a trivial consequence of the 
analysis in this paper. The example of this section has 
been chosen because it pedagogically illustrates well the 
simplest among the ideas of the coming analysis. 



measure, its density r{x) is continuous and < 

fi{S-tdx)/ ii{dx) < for some A > and for all < 
t < 9(x): this is referred to by saying that ji is quasi 
invariant with respect to the stopping time 6{x) on S: 
symbolically ^ is (S, ^^(x))-A-quasi invariant. 

Let dsx be the surface element on S and v{dsx) its 
area computed by measuring the volume with /x. Then 
the following Sinai's lemma, [5, 6, 14], holds: 

Lemma 3: IJ fi is {T,,9{x))-X-quasi invariant the inte- 
gral of any non negative function f over the tube with 
base E and ceiling 6{x) can be hounded by 



[ f{y)lJ-{dy) <e^[ I f{S-rx) \vx ■ n^l v{dSa:)dT, 
Je JuJo 

> e-M / f{S-rx) ■ v{ds^)dT (H.l) 
Jn Jo 



The lemma can been used to reduce dynamical esti- 
mates to equilibrium estimates. 

Proof: Let the trajectory of a point y which reaches E 
within the stopping time at a; e E be parameterized by 

the time r. Then the set of points into which the paral- 
lelepiped A with base dsx and height dr evolves becomes 
a parallelepiped St A. with base Stdsx and the same height 
dr. Therefore the measure of /z(S'tA) is e~^ < < 

hence the integral of any positive function f{y) over 
the set E can be bounded above and below by the inte- 
gral of Jj. J^^^^ f{S-tx)p{x)i'{dsx)dT if p{x)v{dsx)dT is 
the measure of A: the latter is \vx ■ nx\v{dsx)dT. Notice 
that p{x) = \r{x) Vx ■ nx\ if iJt{dx) = r{x)dx. 
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Appendix H: Quasi invariance 



A probability distribution on a piecewise regular 

manifold M is quasi invariant for a flow x — >■ Stx gen- 
erated by a differential equation x = Vx if e-^(*) < 
IJL{S-tdx)/iJL{dx) < e^(*) and \{t) < oo. 

Suppose given > 0, a piecewise smooth surface E C 
M with unit normal vector Ux and a "stopping time" 
X — >■ d{x) < defined on E consider all points a; e E 
which are reached for the first time in positive time t < 
6{x) from data y ^ E. Call E, the set of such points, i.e 
the tube with base E and ceiling 9{x). 

The probability distribution p, is quasi invariant with 
respect to E and to the stopping time x 6{x) if 
it is absolutely continuous with respect to the volume 



Consider N particles in Uj VI j n A with a configuration 
of immobile particles outside A. The analysis in [7] can be 
followed and the solution of the equations of motion can 
be defined on the set r+ consisting of the configurations 
X in which 1 of the particles is at ^ on the boundary 
dK, where elastic collisions take place, with normal speed 
9 ■ "(0 > 0. The time evolution makes sense until the 
time T+{x) of next collision; it can then be continued 
after the elastic collision because, apart from a set of 
zero volume, the normal speed of the collision can be 
assumed 7^ 0, until the time t{^{x) > 0, if any, in which 
the total kinetic energy in one of the containers Ctj, j > 0, 
vanishes. 

Consider the subset r(0) c r+ of the points with 
t{x) < 0. If A C r(0) and if ly is the measure z^o.s'j 
see Eq.(7.9) and appendix D, is a small set it will be 
i/(TA) > A~^i/(A) where A(a;) is an upper bound on the 
entropy production in any interval within [0, r(a;)]. 



J: Appendix: Entropy bounds: check of Eq. (7.7) 



19 



There cannot be infinitely many collisions, except 
on a set oi ly measure, as long as A(T"a;) < L < 
00. If T denotes the return map the only way an 
accumulation of collision times with X^^q t{T'^x) = 

9 < Q could occur is if also K{Sl"''^h) ^0 and 
x)dt = +00 : hence the limit 



supt',t"<e It' '^i^t 
limt^e K{S^"''^^x) = 0. The reason is that if A remains 
finite then Poincare's recurrence argument will apply im- 
plying, as in [7], that J2h=o'''i'^''^) ~ outside a set 
of r^-measure 0. 

This means that the thermostatted time evolution is 
well defined until the first time t\{x) when K {sj.^'^"^ x) 

0. Since the thermostat force 
Kj — )■ the time evolution will be defined in the closed 
interval [0,tA(x)]: the evolution proceeds, well defined, 
until the first time tA{x) (included, if any) when some of 
the K\j vanishes, see remark (6) p. 3. 



tA 

■j^qi is bounded even if 



Appendix J: Entropy bounds: check of Eq.(7.7) 



Writing for the smallest integer > (Xogn)'^ g^{£^/r^p) 
(here is chosen instead of the natural 92^/2 in order to 
simplify the formulae: recall that by defiinition 7 > 1/2), 
then /io.s' almost surely, S' splits into an union over ^ € 
A„ n r^l/ of the union of 5| U 5| U 5^, where 



Si ={y e S' : I2/ n Ql = fcj, \y n dC^\ = 1} 
Si ={?/ e S' : y n Q 9 {q, q), e{q, q)=e^} 

S^ ={y e S' : KJ^r^{y = ^'^2"''} 



(J.l) 



if £j =^ {{lognyg-y{^/r^)). 

Consider first the case of S^. Let D c S^ he the set of 

the X which satisfy Kj^A„{x) = ^k2'^" for a given j > 
while Kj,^AAx) > iK2^""for f > 0, / ^ j. 

Recall the DLR-equations, [26], and consider the clas- 
sical superstability estimate on the cixistence of c > 
such that Pn = e"'^^ " bounds the probability of finding 
more than p2'^'^ particles in A„ n VI j if p is large enough 
{e.g. p > maxj 5j). Then the probability pa{D) can be 
bounded by Pn (summable in n) plus, see Appendix H, 



^^^^ / Mdq') E © ' 
i=i 



(J.2) 



where /Iq is the distribution of the positions in /jq, q = 
{qi, ...,qi)e i^j n A„)', P2 = i«2d», C/A„,,(g, q') is the 
sum of ip{q — q') over pairs of points qi, q'^ € (1 A„ and 
pairs e A„ n flj,q'f ^ A„ fi ilj, and 



(1) ZA„j{q') is the partition function for the region A„ fl 

flj (defined as in Eq.(3.1) with the integral over the g's 
extended to A„ n ilj and with the energies {q, g')); 

(2) The volume element P^'^-^dP has been changed to 



E 



\dq^{q - q' 



E 

qeA„nnj 



(J.3) 



so that PP is a bound on the time derivative 2PP of the 
total kinetic energy contained in A„ evaluated on the 
points of D: hence 2PP is bounded by < CP^P). 

(3) ijj{id) is the surface of the unit ball in M^''. 

(4) The e*-^^^^ takes into account the entropy estimate 
i.e. the bound, CMV , of the non-invariance of /ig. 

The integral can be (trivially) imagined averaged over 
an auxiliary parameter e G [0,e] with e > arbitrary 

(but to be suitably chosen shortly) on which it does not 
depend at first. Then if P is replaced by (1 — £)P in the 

exponential while P*»-^ is replaced by ^ypi^-i-i the 
average over e becomes an upper bound. Changing e to 
Pe {i.e. hence de to (^^^) 2 dPe) the bound becomes the 

;Uo-average ^^^(Px«,?)/i*o bounded by 



25 e 



?C2'' 



£^l/22dn/2 



-c2' 



id/2 



(J.4) 



where Xk.e is tlie characteristic function of the set {(1 — 
e)2|2''" < Kj < f 2'^"}. The inequality is obtained by a 
bound on the first average, via a superstability estimate, 
proportional to 2^''" and by the remark that the second 
average is over a range in which K shows a large devia- 
tion from its average (by a factor 2) hence it is bounded 
above by e"''^"'' with h depending on k but independent 
on £ for n large. Therefore fixing e small enough (as a 
function of n) the bound holds with suitable C, c > and 
is summable in n (and of course on j > 0). 

Similarly, the surface areas i'o,e'(<5|) and i'o,e'(<5|) on 
<S^ , <S^ induced by ^,o are bounded by 



(for suitable C, c, functions of E). Summing (as 7 > 1/2) 
over ^ g j n A„ as discussed below and adding Eq.(J.2) 
the Eq.(7.7) will follow. 

Consider the case of »S|. By Eq.(J.l) if y e »S| then 
|y n CjI = and there is (g, q) G y with q G dC^. 

Remark that y is the configuration reached starting 
from an initial data x G Xe within a time T„(a;) < 9: 
hence Eq.(5.8) applies. By Eq.(5.8) 
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/ Mo,E'(rfy) / 
Jsi Jo 

< eviCn"^' j n{dx) 



dtw{y) 

Ji 
ZcAx) 



(J.6) 



where ^{dx) is the /^o-distribution of configurations x 
outside and 



where the area({i?((7, q) = E{\) is the area of the surface 
{{q, q) : E{q, q) = E^} in R^'^ (the E^ is defined in (J.l)). 
Then J2 is bounded by 

(J.12) 

so that, suitably redefining C, C (functions of E), the 
contribution from <S| is bounded by 



Ji = dqi 



dq2... dqt^ 



dqe 



-l3iH{q,q\x) 



(J.7) 



The estimate of the r.h.s. of Eq.(A.l), as remarked, is an 
"equilibrium estimate" . By superstability, [9] , and since 
<p > 0, the configurational energy U {q\x) > bk^ — b'k^, so 
that Ji is bounded by: 



(fcj-1)! ^Pjm' 



(J.8) 



while J iJ,{dx) Zc {x) — ^ because Zc^{x) > 1: and the 
bound can be summed over k^. Thus the contribution 
from iS| to Eq.(7.7) is bounded by 



(J.9) 



with C, b suitable positive constants. Since 7 > 1/2, this 
is summable over ^ and yields the part of the Eq.(7.7) 
coming from the integration over . 

Let, next, y G iS| and let {q, q) as in (J.l). The function 
w is, if E{q,q) = e{q,qf, 



w 



\dE{q,q)/dt\ 
|grad£(g,g)| 



(J.IO) 



and \dE{q,q)/dt\ < C\q\n^ because dE/dt is the work 
on the particle {q, q) done by the pair interactions (ex- 
cluding the wall forces). It is bounded proportionally to 
the number of particles which can interact with {q,q), 
which, by theorem 4, is bounded proportionally to n^^"^ 
(as the total configuration is in S'). On the other 
hand, |gradi?(g, g)! 



m- 



|gp > m\q\ hence 



w < Cv}/'^ again by Eq.(5.8) and the remark preceding 
Eq.(J.6). 

Then, analogously to (J.6), the integral under consid- 
eration is bounded by Ce^ ^/n {C,C' are suitable 
constants functions of E) times an equilibrium integral 



J ^i{dx) 



J2 

Zc^ (x) 



with J2 defined by: 



(J.ll) 



e-^^^(area{{E{q,q)=E^}) 



C'e 



iCMV 



'ne 2 



ift[(logn)-'g-,(?/r^)]^ 



(J.13) 



and Eq.(7.7) follows from Eq.(J.2), from (J.9) and from 
Eq.(J.13). 



Appendix K: Details the derivation of Eq.(7.5),(7.6) 



To prove item (1) and Eq.(7.6), thus completing the 
proof of theorem 7, we shall compare the evolutions 
with s("'°)(t), at same initial datum x G Xe and 
t < T„(x), the latter being the stopping time defined in 
Eq.(7.4). We start by proving that there is C > so that 
for all n large enough the following holds. 

Lemma 4: For t < Tn{x), see Eq.(7.4), and k > 

(logn)'*', then 



.(n.l) 



\ql--^(t)\<Cvi {k logn)T, 



qr \t)\<r^{2>^ + C{k logn)T). 



(K.l) 



for qi{0) e Afc. 

The necessity of this lemma is due to the fact that we 
cannot control the positions and speed at time t in terms 
of the norms ||x||„ at time 0: since the particles move we 
must follow them (a "Lagrangian" viewpoint). 

A corollary of the above will be: 

Lemma 5: Let J\f and p be the maximal number of par- 
ticles which at any given time < T„(x) interact with a 

particle qi initially in Afe+i and, respectively, the mini- 
mal distance of a particle from the walls. Then 



Af < C (klognf^ , p>c(fclogn) 



-2(d7+l)/a 



(K.2) 



for all integers k G ((logn)''', 2(logn)'''). 



The proof of the lemmas is in Appendix L. The interval 
allowed to k could be much larger, as it appears from the 
proof. 

def 

We have now all the ingredients to bound 6i{t,n) = 

\Qi^'^\t) ~ Qi"'°\t)\- Let fi be the acceleration of the 
particle i due to the other particles and to the walls. 
In the equations of motion for the evolution, Eq.(3.7), 
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for q^"'^\t) the elastic collisions with 9A„ (implicit in 
Eq.(3.7)) can be disregarded if the particle initially in qt 
does not experience a collision with 3A„. For n large 
such a collision cannot take place if qi{0) G ^k+i, k < 
2(logn)'^ and t < T„(a;), by Eq.(K.l). 

Then, for 2(logn)T > fc > (logn)'', by Eq.(K.2) and 
if qi e Afc+i, it follows that \fi\ < C(fclogn)'', with 

rj a suitable constant {ri'^= dj + 2{dj + 1)(1 + ^)), so 
that subtracting the Eq.(3.7) and Eq.(6.1) for the two 
evolutions, it follows that for any qi G Afe+i (possibly 
close to the origin hence very far from the boundary of 
Afe, if n is large) 



number of particles of x "' (t) in A,, which in turn is 
bounded by M up to time O, by theorem 6. An analogous 
argument for the velocities allows to conclude Eq.(7.6). 

Applying again Eq.(7.5) the proof of theorem 7 is com- 
plete apart from the uniqueness statement: the latter 
can be proved (as usual by iteration of the equations in 
integral form) but we skip the details, see [8, 13]. The 
uniqueness is within the solutions of the equations of mo- 
tion in the infinite volume and with data subect to the 
condition of being in Hi/d- 



Si{t,n) < C(fclogn)''2-"<^ 

+ e /* \Mq^-''\T)) - Mq^-''\T))\ dr, ^^'^^ 
Jo 

because, recalling the definition of aj and the third con- 
dition on the stopping time in Eq.(7.2), the denominator 
of |Q!j| is bounded proportionally to 2"''. 
Remarking, by lemma 4 for (/{"'^'(t), that 



Appendix L: Proof of Lemrticis 4,5 

With the notations of Appendix K and 7 > 5: 
Proof (lemma 4): If t < T„(a;) then 

|.(n,i)/+M^ //, X, |ai"'^^(i)l + V2r<^x7 

\ql ' '{t)\ <vi {{log n)log+^ 4^ ^) , (L.l) 



^max^ \qt''\t) - ft| < Cr^{k\ogny, (K.4) 

and, by theorem 5 for q\"'^\t)), also 9-"'°'(t) satisfies an 
identical relation, let f > be an integer, ke such that 



2fe, =2'=+^C(fclogn)^ 



(K.5) 



and let Uki{t,Ti) the max of Si{t,n) over \qi\ < r^2''^. 
Then by Eq.(K.3) and Eq.(K.2), 



Uke{ t,n 
r, 



^ < C(fclogn)''(2-"'^+ / 
V -'0 



* Uki^^{s,n)ds 



(K.6) 

for i<t= 2'=/((fc log n)''C), the latter being the largest 
f such that 2'=^ < 2'=+!. By Eq.(K.6),(6.10), and Eq.(K.2) 

^ ^cik ^o^nT c{k\ogn)^2-''^ (K.7) 

for (logn)^ < k < 2(logn)'''. 

Thus Uk{t,n) is bounded by the r.h.s. of the first of 
Eq.(7.5); analogous argument shows that also the veloc- 
ity diff'erences arc bounded as in Eq.(7.5) which is thus 
proved for all t < Tn{x). 

Therefore given qi{0) with \qi{0)\/rip < 2*"' it is, for n > 



e'^o^^ large enough and i fixed, \ql"'''^\t) — qi"''^\t)\/r^ < 
,)-,{t,n) < C2-^(i°s")\ i.e. for n large '^^(i) is 



i(log 



closer than to q^/^'^\t). 

It remains to check Eq.(7.6): for n large the number of 

particles in x\"'^\t) which are in A* is smaller than the 



implying: |9|"''^^(i)| < r{t)r^ if r(f) is an upper bound 
to a solution of Eq.(L.l) with = replacing < and ini- 
tial datum \q^'^'^\0)\ < 2'=r^. And r{t) can be taken 
r{t) =^ 2'= -I- 2t;i ((log n) log+ 2*=)^ for t < 6, provided 

((logn)(log+r(e) + V2)y < 2((logn) log+ 2'=)^ (L.2) 

which, since (*^^°s") vanishes as n diverges (because 
k > (logn)''), is verified for all n large enough. Thus 
\q'r'^\t)\ < r^r{t), hence \q'r'^\t)\ < r^Cf{t) for all 
t < Tn{x), and the lemma is proved. 

Proof (lemma 5): By lemma 1 the following properties 
hold for all n large enough and all t < Tn{x): 

(i) for all qi G Afe+2, 



max \qi"'^\t) - qi\ < C r^{klogny , 



(L.3) 



which implies the first of Eq.(K.2). 

(ii) particles in A^ do not interact with those 

By Eq.(L.3) we see that if qi e Afe+i then q^i^'^\t) € 
Afc+2 so that, by the definition of the set B, recalling that 
7 > 1/2, the lemma follows because the work of the pair 
forces on qi is bounded A/" < C(fc logn)''''' times the speed 
bound C(A:logn)'', by Eq.(K.l) and the second of (K.2). 
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